Signatures of sub-GeV dark matter beams at neutrino experiments 



Patrick deNivervilleQ David McKeenQ and Adam Rit^ 
Department of Physics and Astronomy, University of Victoria, Victoria, BC V8P 5C2, Canada 

(Dated: March 4, 2013) 

We study the high-luminosity fixed-target neutrino experiments at MiniBooNE, MINOS and T2K 
and analyze their sensitivity to light stable states, focusing on MeV-GeV scale dark matter. Thermal 
relic dark matter scenarios in the sub-GeV mass range require the presence of light mediators, 
whose coupling to the Standard Model facilitates annihilation in the early universe and allows for 
the correct thermal relic abundance. The mediators in turn provide a production channel for dark 
matter at colliders or fixed targets, and as a consequence the neutrino beams generated at fixed 
targets may contain an additional beam of light dark matter. The signatures of this beam include 
elastic scattering off electrons or nucleons in the (near-) detector, which closely mimics the neutral 
current scattering of neutrinos. We determine the event rate at modern fixed target facilities and 
the ensuing sensitivity to sub-GeV dark matter. 



1. INTRODUCTION 

The existing gravitational evidence for dark mat- 
ter provides limited information about its non- 
gravitational interactions, and many candidates are 
sufficiently non-relativistic and weakly interacting. 
The paradigm of a weak-scale thermal relic has the 
virtue of simplicity, with an abundance fixed without 
detailed knowledge of early- universe physics. How- 
ever, direct detection experiments now impose strin- 
gent constraints on dark matter with a weak-scale 
mass; for example, spin-independent cross sections 
on nucleons must be at or below 10~^^ cm^. With 
this sensitivity now crossing the Higgs-mediation 
threshold, the minimal weakly-interacting massive 
particle (WIMP) paradigm may need generalization 
to allow new interaction channels, beyond the elec- 
troweak sector of the Standard Model (SM). This 
would position dark matter as part of a more com- 
plex hidden sector containing additional light states. 
The required relic density could then be achieved 
without either weak-scale interactions or a weak- 
scale mass [iHi]- 

This viewpoint has some interesting implica- 
tions when one looks at the existing limits on di- 
rect WIMP scattering. The sensitivity of direct- 
detection experiments tends to fall rather sharply 
for masses below a few GeV, due to the recoil en- 
ergy detection threshold. The GeV mass scale also 
happens to coincide with the Lee- Weinberg bound 
[5], below which a thermal relic needs non-SM anni- 
hilation channels through light states to ensure the 
correct relic abundance. In combination, these ob- 
servations naturally lead us to explore the use of 
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new experimental tools to probe the sub- GeV mass 
range for thermal relic dark matter. The presence 
of light mediators coupled to the SM opens up the 
possibility of producing these states directly in ac- 
celerators or fixed target facilities. This 'dark force' 
phenomenology has been the focus of considerable 
interest in recent years. For example, a number of 
search strategies are based on the production of a 
GeV-scale vector mediator, with its subsequent de- 
cay to lepton pairs [2^ l6Hllj. However, these search 
strategies are limited if, instead, the mediator is not 
the lightest hidden sector state and decays predom- 
inantly into the hidden sector, e.g. to dark matter. 
In this case, the scattering of those light states in 
a detector spatially separated from the production 
point represents perhaps the most efficient search 
strategy. Moreover, owing to the potentially large 
production rate, and the existence of large volume 
(near-) detectors, proton fixed-target facilities focus- 
ing on neutrino physics appear to be an ideal means 
for exploring these scenarios. 

In this paper, we analyze the sensitivity of neu- 
trino facilities to a boosted light dark matter beam 
produced via the generation and subsequent decay of 
GeV-scale mediators. This extends our earlier anal- 
ysis of MeV-scale dark matter [11] [12] to the full 
sub- GeV range. We will find that high-luminosity 
experiments such as MiniBooNE, MINOS and T2K 
have significant sensitivity to neutral current-like 
scattering of sub-GeV dark matter off nuclei in the 
(near-) detector. Although there is a long history of 
searches for exotics using fixed target facilities (see 
e.g. pi [BHTTj [T3] ) , neutrino experiments have the ad- 
vantage that the large detector volume is sensitive 
to scattering signatures in addition to the products 
of SM decays. Since the recoil energy of sub- GeV 
halo dark matter is generally below threshold for un- 
derground direct detection experiments, and search 
channels at high energy colliders are less sensitive 
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in the case of light mediators, high-luminosity fixed- 
target experiments can play a complementary role 
in direct searches for dark matter. 

In order to be as model-independent as possi- 
ble, we parametrize the mediator interactions via 
the lowest dimension operators (portals) for a SM- 
neutral hidden sector, Ci^t = "^OsmOhs^ where O 
denotes SM and hidden sector (HS) operators. For 
light dark matter, fixed-target facilities have an ad- 
vantage if the mediator can be produced on-shell, so 
we focus on the renormalizable vector [I4 and scalar 
SM portals [15]: 



Ant = >Chid(-^, X) 



hiF^j^V^^j vector portal 
ASH^H^ scalar portal 



(1) 

where and H are the hypercharge field strength 
and the Higgs doublet, while >Chid provides hidden 
sector couplings between the mediator field X = 
or S and the light dark matter candidate x- We will 
limit attention to the kinematic regime 



mx > 2m^ - 0{MeV - GeV), 



(2) 



so that with small portal couplings to the SM, the 
mediators predominantly decay into the hidden sec- 
tor, Br(X ^ xx) ^ 1- 

The rest of this paper is organized as follows. In 
Sec. [2j we describe our model for production of the 
dark matter beam at MINOS, T2K and MiniBooNE 
using both vector and scalar portals. In Sec. |3j we 
discuss a number of existing constraints on sub-GeV 
dark matter, coupled to the SM via these portals, 
detail the annihilation and scattering rates, and de- 
termine viable models which can be probed using 
neutrino facilities. In Sec. [4j we focus on the most 
viable dark matter scenario, with scalar dark matter 
coupled via the vector portal, and analyze the sensi- 
tivity to the ensuing dark matter beam at MINOS, 
T2K and MiniBooNE. We conclude in Sec. O 



2. PRODUCTION OF THE DARK MATTER 
BEAM 

A. DM interactions 

The viability of thermal relic dark matter with a 
mass in the MeV-GeV range, well below the Lee- 
Weinberg bound, seemingly rests on the presence of 
a light hidden sector with states that can mediate 
annihilation [2-[4]. Moreover, various phenomeno- 
logical constraints [4] suggest that the most viable 
scenarios are those in which the hidden sector is un- 
charged under Standard Model symmetries. This 
naturally leads us to the portal interactions ([T]) as 



the primary means of probing these sectors at low 
energies. 

To keep our analysis as general as possible, we 
will consider both the vector and Higgs portals for 
production of the dark matter beam in this section. 
These light mediators are necessary to allow for a 
viable annihilation channel in the early universe, but 
we will be agnostic about the precise choice of model. 
This will allow us to analyze the raw sensitivity of 
neutrino facilities to production of these light states, 
and we will turn to the model-dependent constraints 
on viable light dark matter scenarios in the following 
sections. 

To fix the interactions, we use the simplest realiza- 
tions for the vector and scalar portals, and moreover 
we will only need their low energy manifestations. 
For the vector portal coupling, F^V^^, we have 



Cv = {eU^^ + eV^) + C^^^iV, x) 



(3) 



where we have used d^F^^ = eJ^^ in terms of the 



■ ■ , while 



electromagnetic current J^^ 

is the corresponding U(l) current for scalar (or Dirac 
fermion) dark matter, with gauge coupling e'. 



jt^ = I ^X'^^X + 0{V^), scalar 
X \ ^XT^X? fermion 



(4) 



^kin{V^ x) contains canonical kinetic and mass terms 
for V and x^ ^ind higher order potential terms have 
not been written explicitly. 

For the trilinear scalar portal coupling, SH^H^ we 
have 

£3 = 3 (^JgWsB + 13 J;^) + Ain(5, X) + • • • (5) 

where we have integrated out the SM Higgs, which 
induces a coupling ~ Av/mj^ between S and the 
SM fermions via Jewsb = ^qQQ/^ ~^ ' ' ' ^ while 
is an analogous mass current for scalar (or fermion) 
dark matter. 



jm _ J ^xX'^X, scalar 



^XX, 



fermion 



(6) 



We have inserted a factor of in the scalar case, so 
that 13 remains a dimensionless coupling. As above, 
^kin{V^ x) contains the kinetic and mass terms for S 
and scalar /fermion dark matter x- 

We will refer to the mediator V or S generically as 
X, and the crucial kinematic assumption will be that 
mx > 2m^, so that the mediator can decay on-shell 
to dark matter. For small mixing via the portals, 
the hidden sector branching Br(X xx) ^ 1- 



B. Production mechanisms 

There are two viable production mechanisms for 
the mediator X at proton fixed-target experiments: 
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• Direct production: This corresponds to 
hadron-level processes such as pp{n) X* 
XX (or X^x) shown in Figs. 1 and 2. In prac- 
tice, since X can decay to xx, we will use the 
narrow width approximation so that X is pro- 
duced on-shell. In this approximation, valid to 
O (e^^, , the cross section for the production 



of a DM pair can be written as 

cr {pp{n) ^ X* ^ xx) 

= a{pp{n)^X)BT{X^XX)- (7) 

The direct production cross section of a vector 
mediator is 



a {pp{n) V) 



i: 



dx 



do- {pp{n) V) 



dx 



dx 



fq/p {X) fq/p^ri) (^) + fq/p (^) fq/pM (^) 



(8) 



where Cq is the charge of quark q in units 
of the positron electric charge, r = rriy/s, 
and ^/s is the hadron-level center-of-mass en- 
ergy. The parton distribution function (PDF) 
fq/p{n) (^) gives the probability of extracting 
the quark q with momentum fraction x from a 
proton (neutron) and similarly for fq/p(^n) (^)- 
We have omitted the scale, Q, at which the 
PDFs are evaluated. To obtain estimates, we 
use CTEQ6.6 PDFs [16] and set Q = my; 
varying Q in between my /2 and 2my resulted 
in an uncertainty in the production cross sec- 
tion of less than 30% for my > 1 GeV at 
T2K and MINOS beam energies. Higher-order 
QCD corrections are large, introducing an er- 
ror that can potentially be 0{1). 

The production cross section as a function of 
the DM lab frame energy, E^^ and the an- 
gle between its lab frame momentum and the 
beam direction, 0, can be related to the differ- 
ential cross section in Eq. (Isl) through 



da {pp{n) ^ V ^ xx) 
dE^d cos 



9(x, cos 6) 
d{E^,cose) 



(9) 



^ daippin)^V) B,(^^^^),(eos^), 

where is the angle between the momentum 
of X and the beam in the V rest frame and 
the quantity in square brackets is the Jaco- 
bian associated with this variable change. The 
function g describes the angular distribution 
of the DM in the V rest frame. For scalar DM 
produced through a vector mediator, this is 



g (^cosO^ = ^ (l-cos^O^ . 



(10) 



If, instead, x is a Dirac fermion, then 

g(^cose^ = ^ (l^cos^O^ . (11) 



We will find the distribution of V momenta 
useful. 



fv {pv) 



1 



da {pp{n) V) 



a {pp{n) V) dpy 

1 dx da {ppin) V) 



(12) 



a {pp{n) V) dp\ 



dx 



with py the momentum of V in the lab frame 
which is related to x through 



Pv 




(13) 



1/2- 



{x — r/x) p'^mj 



where mr = ^p,n is the target mass, is the 
momentum of the beam, 7/3 = Pb/Vs, and 

7 = 

For illustration, in Figs. |3]-[5] we present the 
resulting direct production distributions for a 
vector mediator that subsequently decays to 
scalar DM at the T2K and MINOS experi- 
ments, where £^beam = 30, 120 GeV {^/s 
7.6, 15.1 GeV), respectively; see Sec. [4| for fur- 
ther details of these experiments. Fig^ shows 
the total production cross section for pp and 
pn collisions at T2K and MINOS as a func- 
tion of the vector mediator mass. After inte- 
grating over energy, the angular distribution 
of scalar DM is shown in the top of Fig. [4] 
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FIG. 1. Direct production of scalar dark matter via the vector portal. The leading-order process is shown on the left, which 
is helicity suppressed in the forward direction. The process on the right is higher order in a^, and also phase space suppressed, 
but has less helicity suppression in the forward direction. 




FIG. 2. Direct production of dark matter via the scalar por- 
tal. The solid gluon fusion ggh vertex is generated at 1-loop. 
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FIG. 3. The total production cross section of a vector media- 
tor at T2K and MINOS energies as a function of the mediator 
mass. The solid and dashed curves indicate the cross sections 
for pp ^ V and pn ^ V respectively. The PDF scale has 
been fixed to Q = my- 



at T2K ND280 and MINOS in the case that 
mv = 1 GeV and = 300 MeV. We fo- 
cus on the off- axis ND280 detector at T2K, 
to contrast with the on-axis detector at MI- 
NOS in sampling the angular production dis- 
tribution. However, comparing ND280 to the 
on-axis INGRID detector at T2K would pro- 
vide a similar contrast. In the bottom left of 
Fig. |4j we zoom in on the relevant angular re- 
gion for the off-axis T2K ND280 near detector 



and show the scalar DM angular distribution 
for mv = 1 GeV and several DM masses pro- 
duced in pp collisions. We do the same in the 
range of angles around the MINOS near detec- 
tor in the bottom right of Fig. [4j As the mass 
of the DM is increased, it is produced in the 
more forward direction since its velocity in the 
V rest frame decreases. However, the angu- 
lar distribution of scalar DM produced via a 
vector mediator, Eq. (10), suppresses the pro- 



duction of DM along the beam direction itself. 
Thus, despite the smaller cross section for the 
production of vector mediators as a result of 
the lower energy of its beam, a larger number 
of DM particles may pass through the off-axis 
T2K ND280 near detector than the on-axis 
MINOS near detector. This suppression along 
the beam axis is lessened somewhat when con- 
sidering higher-order production mechanisms 
like the diagram on the right of Fig. [T] which 
we do not include in this study. 

We show the energy distribution of scalar DM 
for mv = 1 GeV and a range of m^ in pp 
collisions for T2K at l9 = 2° and for MINOS at 
= 0.025° in Fig. ^ 

For a scalar mediator, the leading-order direct 
production cross section is 



a {pp{n) S) 



2/)2 



288A/27r 



(14) 



Here, r = m^/s and the PDF fg{x) is the 
probability of finding a gluon with momentum 
fraction x in a nucleon. Up to threshold effect 
corrections, N counts the number of quarks 
with a mass greater than ~ 0.2m5' [17 . 

The DM distributions in the lab frame can be 
related to the differential production cross sec- 
tion in the same way as in the vector mediator 
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FIG. 4. Top: production distributions of scalar DM as a function of lab frame angle with respect to the beam, normalized 
to unity, in the case of a vector mediator with vfiy = 1 GeV. The solid curves indicate pp collisions and the dashed ones 
pn collisions. The set of curves that peak at ^ ~ 1° corresponds to a p beam with an energy equal to that of the MINOS 
experiment (£^beam = 120 GeV, = 15 GeV) while the curves that peak at ^ ~ 4° correspond to a beam with an energy 
equal to that of T2K ND280 (£^beam = 30 GeV, = 7.6 GeV). Bottom left: production cross sections in the case of a vector 
mediator (mv = 1 GeV) and scalar DM for m^ = 100, 300, 450 MeV (solid, dashed, dotted) as functions of the DM angle 
with the beam in the lab frame in the case of pp collisions at an energy corresponding to the T2K experiment. The range of 
angles shown coincides with those covered by the off-axis ND280 near detector at T2K. Bottom right: the same at MINOS 
beam energy. The angles shown here are those that the MINOS near detector covers. 



case in Eq. (|9|. In the S rest frame, the DM 
is simply produced isotropicahy, 



g{cosO) 



1 



(15) 



Because of the weak scale and loop factor sup- 
pressions in Eq. (14), scalar mediator produc- 



tion is extremely small compared to the that of 
a vector at GeV scales. Thus, current neutrino 
experiments are much less sensitive to DM sce- 
narios involving a GeV-scale spin-0 mediator 
than they are to a spin-1 mediator. Along with 
other factors to be discussed in the next sec- 
tion, this will lead us to focus only on the direct 
production of vector mediators. 



• Indirect production: This corresponds to pro- 
duction of X via the decay of hadronic states 
(generically denoted (j)) produced in the pri- 



mary pp and pn interactions, 
p^p(n) ^ 



X 



XX 



(16) 



Depending on the beam energy and form of the 
target, the relevant decay lengths ensure that 
this entire sequence of events will occur either 
inside the target itself or in the subsequent de- 
cay volume. 

In practice, this process is most important in 
the low mass range where, for example, the 
large production rate of neutral pseudoscalar 
mesons (j) = 7r^,?7 can dominate the over- 
all production of F's in particular. The me- 
son production distribution at MiniBooNE is 
well-described by the Sanford and Wang fit 
f^^{0,p) as described in [19 , and utilized pre- 
viously in [12^. To estimate the meson pro- 
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FIG. 5. Left: da {pp ^ V ^ XX^) /dE^dO in the case of scalar x ^-nd vector V for mv = 1 GeV and = 100, 300, 450 MeV 



(solid, dashed, dotted) for pp colhsions at T2K ND280 at 



2°. The cusps at the kinematic limits for larger are the 



result of a degeneracy in the angle between x ^i-nd the beam direction in the lab frame, ^, as a function of its value in the V 
rest frame, ^, for relatively small DM velocities (e.g. in the limit that the DM is produced at threshold, ^ = for all 6). Right: 
The same for pp collisions at MINOS at ^ = 0.025°. 
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FIG. 6. Pseudoscalar meson production distributions f^^^^ {0,p) in angle (left) and momentum (right), according to the 
fit J8], scaled to the beam energy and target composition for MINOS and T2K. These distributions determine the indirect 
production of vectors via e.g. rj — V^. Note that the fitted distributions have an unphysical low-energy tail at large angles 
which we exclude in the analysis. 



duction distribution for MINOS and T2K, we 
make use of an analytic fit f^^^^{0,p) [H] to 
data for (averaged 7r+ and 7r~) pion production 
obtained over a range of energies, which can be 
scaled to cover the target materials for the ex- 
periments of interest. Example distributions 
for MINOS and T2K configurations are shown 
in Fig. |6] We have also tested this distribution 
against existing data published by NA61 [20^ 
for the T2K target configuration, and found 
good agreement. The aforementioned NA61 



data is used to estimate the total pion flux at 
T2K. As there is currently no equivalent pion 
production dataset for MINOS, the pion pro- 
duction cross section was estimated by scaling 
the NA61 measured total cross section using 
the relative magnitudes of the BMPT distri- 
butions for T2K and MINOS Q The produc- 
tion rates of negatively and positively charged 



The total 7r+ fiux could also be determined by working 
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pions are averaged, and differ by 0{1) fac- 
tors (see e.g. [18 ). To estimate the r] pro- 
duction rate, we use this averaged distribution 
and make use of some early experimental data 
[2T| [22] , which indicates that in the appropri- 
ate energy range 

^pp^ppn (25 30) CFpp^ppr]' (-^'^) 

This production mode is most relevant for V, 
and in the case = tt^ or 77 the branching ratio 
to V is proportional to that of the radiative 
decay to two photons, though suppressed by 
coupling and phase space factors, 

Br^^^y ^ - ^ j Br^^^^. (18) 

For 77 decays, as will be relevant here, 
Br^^^^ 0.39. This process becomes less 
competitive for higher mediator masses. We 
also explored the rate of V production due to 
radiative decays of cc mesons such as J/^l) 
Vrj, but the overall rate is well below that of 
direct production discussed above. 

After either direct or indirect production, the sup- 
pressed portal couplings ensure that the real or vir- 
tual X has an order-one branching to the hidden 
sector, Brx^2x — 1- For completeness, we note that 
away from thresholds, the ratio of the decay rate of 
V to a single SM state [10] relative to the hidden 
sector is of order K^aja' ^ but is enhanced near reso- 
nances in hadronic channels and by the larger num- 
ber of final states. In practice, Br(V XX) = 1 to a 
good approximation, apart from a small region near 
threshold of size (1 — 4m^/my)"^ < n^a/a' where 
the exponent n = 1/2 (fermionic x) n = 3/2 
(scalar x)- We will excise this near-threshold region 
from the mass range. 

The dark matter beam then propagates along with 
the neutrinos. For the couplings considered here, it 
has a weak-scale scattering rate with normal matter 
and is detectable through neutral current-like elastic 
scattering processes with electrons or, of most rele- 
vance here, with nucleons. We will utilize the param- 
eters and datasets of MiniBooNE, MINOS and T2K 
to probe this scenario. Importantly, MiniBooNE 
has published an analysis on neutrino elastic scatter- 
ing, which DM scattering will closely mimic, which 



backwards from the measured neutrino flux, accounting for 
the angular acceptance of the detector [J. However, this 
reconstruction is more complex when the majority of pions 
decay in flight and are affected by magnetic focusing horns. 



allows some estimate of backgrounds and efficien- 
cies. MINOS and T2K allow access to a higher mass 
range. We employed a simulation to determine the 
dark matter flux incident on the detector, which will 
be described in more detail in Sec. |4j after we have 
considered the viable model scenarios. 



3. MODEL SCENARIOS 

The preceding analysis is applicable to generic sce- 
narios of hidden sector states, coupled through the 
vector and scalar portals. In this section, we will 
study more concrete models of sub-GeV dark mat- 
ter, where the mediator mass satisfies mx > 2m^, 
focusing on thermal relics for which the abundance 
provides a constraint on the annihilation cross sec- 
tion. 



A. Constraints 

The constraints that we will take into account in- 
clude: 

• Relic abundance: If x is a thermal relic, 
the WMAP constraint on the relic density 
^DM^^ ^ 0-1 ^ (0-1 pb)/(cri')fo constrains 
the annihilation cross section at freeze-out to 
be {crv)fo ^ Ipb. Even if the hidden sector 
state does not provide the dominant contribu- 
tion to dark matter, the overclosure constraint 
{o'v)fo > Ipb applies more generally. 

• Impact on the CMB: In order not to dis- 
tort the CMB due to energy injection into 
the IGM through annihilation, there are again 
restrictions on the annihilation cross section 
that become particularly severe for light DM. 
The constraint takes the form f{z){av)cMB < 
0.1 (mDivi/GeV) pb, where f{z) is a redshift- 
dependent efficiency factor which for low mass 
varies from / ~ 0.2 for pion to / ~ 1 for 
electron final states [23^. These limits essen- 
tially exclude a thermal relic below a few GeV 
with an abundance fixed via 5-channel annihi- 
lation. Even for asymmetric DM, the require- 
ment that the symmetric component annihi- 
late away efficiently prior to decoupling leads 
to a lower bound on the annihilation rate of 
the same order [24]. This narrows down the 
field of viable models to those with (velocity- 
suppressed) annihilation since v ^ 10~^ in this 
epoch. Other indirect signatures of annihila- 
tion in the galaxy (where v ^ 10~^) are then 
necessarily suppressed as well. 
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• Visible decays: Focusing first on the vector 
portal, we note that models where V decays 
predominantly to the dark sector are less con- 
strained than those in which V is metastable 
and decays mainly to the SM. For example, 
the fixed target constraints on dark forces via 
leptonic V decays |8l lllj are avoided here for 
this reason as V decays promptly to the dark 
sector. However, there are constraints from 
high-luminosity colliders (particularly the B- 
factories in the case of GeV-scale vectors) 
which are sensitive to rare, but prompt, V de- 
cays to the SM. In comparison to dark force 
searches where in the appropriate mass range 
Br(V ^ - this SM branch- 
ing is suppressed here by an additional fac- 
tor of 0{ahz'^ /a^). A dedicated analysis for the 
higgs'strahlung signature [10 was recently car- 
ried out at BaBar [25 , leading to limits trans- 
lating to K,^ < few X 10~^ in the present sce- 
nario. However, this only applies when the 
dark Higgs is heavy enough to decay predom- 
inantly to two Vs. Although there are no 
specific analyses, one can infer 5-factory (and 
^-factory) limits on generic continuum pro- 
cesses, e+e~ V*7 (and any ex- 
clusive decays, e.g of T{nS) or (/), not forbid- 
den by C-parity [26]). The lack of significant 
peaks in similar rare-decay analyses suggests 
that limits, which in the present case trans- 
late to {a/a')f<i^ < 10~^, apply more generally 
019]. 

For the Higgs portal, there are significant B- 
factory limits on rare B decays, Bt{B K -\- 
^) < 10~^, which directly constrain DM cou- 
pled via the Higgs portal [27| due to decays of 
the form B ^ K -\- x -\- X- 

• Invisible decays: The scenarios considered here 
allow for the mediator to decay on-shell to dark 
matter, leading to new invisible decay channels 
[2] |28l . Searches for rare radiative decays 
with missing energy, such as J/?/^ ^ 7 +^ [30] , 
are limited in the present case by C-parity but 
there are also generic limits on purely invis- 
ible decays of J/ip and T(15') [31 that con- 
strain e.g. J/ip ^ V* ^ XX^ ' Off-resonance, 
the limit is of order a'hi'^ < few x 10~^ [29] , 
which is weaker than the visible decay con- 
straint for a' ^ a but becomes more signifi- 
cant for larger values of a' ~ 1. Moving close 
to the resonance, where e.g. my ^ ^j/^, the 
limit becomes particularly stringent, jo' < 
few X 10~^. However, since the V is still quite 
narrow for perturbative values of a', this only 
applies in a small V mass range, which for J/tJj 



and T(IS') decays is somewhat above the scale 
considered here. We note that future limits on 
invisible decays, e.g. of 0, would be sensitive 
to these scenarios. 

• Energy injection during BEN: Energy injec- 
tion from decays of GeV-scale mediators in 
the early universe can be problematic if it oc- 
curs through hadronic channels before about 
10~^ s. In the present case, most decays will 
occur to the hidden sector. Furthermore, the 
couplings are sufficiently large to ensure that 
decays occur well before BBN, so there are no 
significant constraints from this source. 

• Self -interactions: The models we are consider- 
ing here will become nonperturbative when the 
couplings exceed the naive dimensional analy- 
sis (NDA) scale of order a' ~ 47r (for the vec- 
tor portal) and j3 ^ 1 (for the scalar portal). 
These rough limits characterize the point be- 
yond which our perturbative analysis breaks 
down, with a transfer cross section of order 
c^trans ~ ^iT [a' ^ TTi^ / m\ ^ 10~^^cm^ for a 
GeV-scale mediator. However, physical con- 
straints on self-interaction are comparable in 
the low mass range. Limits on halo ellipticity 
generally require that the transfer cross section 
for scattering is below about 10~^^ — 10~^^ cm^ 
[32] : some limits down to 10~^^cm^ do appear 
in the literature (see [24l [J3 for recent discus- 
sions). This leads to similar limits on a' and 
13 as the NDA limits quoted above. 

B. Models and annihilation rates 

In this subsection we briefly outline the param- 
eter space for a couple of DM scenarios, with the 
above constraints in mind. We will focus on thermal 
WIMPs, so that the annihilation rate determines the 
relic density. With sub-GeV masses, there are strin- 
gent constraints on models with 5- wave annihilation 
cross sections, so we will focus on the cases that 
are p-wave suppressed, namely scalar DM coupled 
to the vector portal and Major ana DM coupled to 
the scalar portal. 

• Scalar dark matter with a U(l) mediator. The 
model contains four parameters; the masses 

and my of the dark matter candidate and 
the vector mediator, the U(l)' gauge coupling 
e', and the kinetic mixing coefficient On re- 
quiring that X comprises the majority of dark 
matter, the constraint on its relic abundance 
allows us to fix one relation between these 
four parameters. The primary quantity here 
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is the annihilation rate, which for the GeV 
mass range of interest is given by s-channel di- 
agrams with e+e~, and hght hadronic 
final states. In the limit of small mixing, we 
can approximate this rate by 

(cr^)ann,V ^ (crv)e + {(^v) ^(l + R{s = 4m^)), (19) 

where R = cr e+ e- ^hadrons / e+ e- ^ fi- , and 
the leptonic annihilation rate is [3 , 



\^^)/ o / 



(20) 



2mi 



m 



2 • 



For <C my ~ 1 GeV, this rate scales 
as {gv)i - 10-33/^2(a7a)(m^/100MeV)2cm2 
using ~ 0.3 at freeze-out. Accounting for all 
annihilation channels, the observed relic den- 
sity r^DM^^ ^ 0.1 ~ (0.1 pb)/(cri')fo reduces 
the number of free parameters to three via a 
constraint that fixes a' = a' {m^^mv ^ i^)- 

The p-wave suppression of annihilation for low 
velocities allows this process to satisfy the 
CMB constraints [2F alluded to above, as well 
as the galactic annihilation fiux limits. 



• Majorana dark matter with a scalar media- 
tor: This model is a natural hidden sector 
generalization of the minimal model of scalar 
dark matter, and there are again four pa- 
rameters: the masses and m^, the hid- 
den sector coupling P and the mixing angle 
0. The abundance constraint will again al- 
low us to determine e.g. /3 = ^{m^^ms^O). 
Annihilation proceeds in the p-wave via mix- 
ing with the Higgs in the 5-channel, and thus 
the rate is dictated by the light Higgs width. 
Given that ms <C m/^, the cross section scales 

as (crv)ann,s ^ p'^ O'^ m\{v'^) {T / TTih*) / {ml - 
4m^)^, where /i* refers to a virtual Higgs of 
mass 2m^. This rate is Yukawa- suppressed 
by the Higgs width F^* unless the Higgs mix- 
ing angle ~ Av/m\ and {3 are relatively 
large. However, these couplings are in turn 
constrained by the 5-factory limits on rare B- 
decays with missing energy. The limits of [27] 
imply that /3^6>^m^v|^/m| < 0(1), which 
we see is quite stringent for mediator and DM 
masses in the GeV range. These 5-decay lim- 
its make this scenario quite problematic as a 
model of thermal relic dark matter. 



C. Summary 

We can summarize the conclusions as follows for 
the four scenarios covered here: 

• Vector portal, scalar x- This DM candidate 
exhibits p-wave annihilation, which is crucial 
to satisfy galactic and CMB annihilation lim- 
its, and thus is viable for sub-percent mixing 
via the portal coupling. The sensitivity of neu- 
trino facilities is significant in this case. 

• Vector portal, fermionic x- This implies 5- 
wave annihilation, and the CMB annihilation 
limits can only be satisfied if x is a highly sub- 
dominant component of WIMP dark matter, 
or has a more complex thermal history. 

• Scalar portal, scalar x- This again implies 5- 
wave annihilation, suppressed in this case by 
the small Yukawa couplings. Thermal freeze- 
out would necessitate large mixing, which 
would rule out this scenario either due to the 
limits on rare 5-decays, or its impact on the 
CMB. 

• Scalar portal, fermionic x' This DM candi- 
date exhibits p-wave annihilation, which can 
avoid the CMB constraints, but to ensure the 
correct relic abundance the mixing must again 
be large, which is strongly constrained by rare 
5-decays. 

Combining this information with the knowledge 
that DM coupled via the Higgs portal has a sup- 
pressed production rate, we conclude that rare B- 
decays provide a more sensitive probe of the Higgs 
portal than dark matter beams. However, scalar DM 
coupled via the vector portal is a viable model and 
we will focus on this scenario in the next section, 
where we outline the sensitivity of neutrino facilities 
to a GeV-scale dark matter beam. 



4. SENSITIVITY TO A GEV DARK 
MATTER BEAM 

A. Scattering 

The detection strategy studied here uses elastic 
scattering of the DM beam in the (near-)detector. 
We outline below the relevant cross sections, focus- 
ing on the vector mediator for the reasons discussed 
in the previous section. 

• Vector-mediated scattering: For a vector medi- 
ator, the scattering of scalar DM on nucleons 
shown in Fig. [7| is similar to neutrino- nucleon 
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elastic scattering (see e.g. [S]), and the cross 



section takes the form 



■XN 



a'K^ 47ra^[Flj,{Q^)A{E,E^ 



\Fl^{Q^)B{E,E^)\ 



dE^ 



{ml + 2mN{E - E^)y {E^ 



(21) 



where E and E^ are the energies of the inci- 
dent and outgoing dark matter particles, re- 
spectively and = 2mN{E — E^) is the mo- 
mentum transfer. We use simple monopole 
and dipole form-factors, Fi^n = Qn/{^ + 
QVm^)2 and F2,n = ^n/{1 + QVrn%)^, 
where Qp = 1^ Qn = 0^ Hip = 1.79 and /^^ = 
— 1.9. The functions A and B are defined as 



A{E,E^) = 2mNEE^ 
B{E,E^) = {E^-E) [{E 



ml{E-E^), 



■Ef 



(22) 
(23) 



-V2mN{E^-E)-Am\] 



B. Simulation 

A Monte Carlo simulation was employed to deter- 
mine the kinematics of both the directly and indi- 
rectly produced dark matter beams. With the kine- 
matics in hand, it is possible to calculate the ex- 
pected sensitivity of MiniBooNE, MINOS and T2K 
to the hidden sector scenarios discussed in Sec. [Sj 
and the production channels of Sec. [2] A detailed de- 
scription of the simulation of the indirect production 
channel at MiniBooNE can be found in [12 . Some 
of the pertinent parameters for each experiment are 
listed in Table [T| and we include some additional re- 
marks below: 



• Scalar-mediated scattering-. For a scalar me- 
diator, the t-channel scattering cross section 
of Majorana fermion DM on nucleons takes a 
similar form to the monopole contribution to 
(21), but in place of H?aa' the cross section 



is suppressed by a factor /^^^^^at/ts/'^eW' 
where we have dropped isospin-violating cor- 
rections to the Higgs-nucleon coupling, and 
retained just the dominant contribution from 
fTs - 0.118 where m^/r, = {N\mqqq\N). Up 
to scalar mixing, this is analogous to conven- 
tional Higgs-mediated scattering and thus is 
quite suppressed relative the vector case above. 
Given the suppressed production rate, we will 
not consider this case further in this section. 



To account for the isotopic content of the detector 
material, we will use an effective differential cross 
section, given by 



Z day 



A - Z day 



dE 



f 



A dEy 



dEy 



(24) 



This expression is an approximation which ignores 
the differing detection efficiencies for scattering of 
bound nucleons. However, the Q^-dependent effi- 
ciency factors quoted by MiniBooNE [19 , for exam- 
ple, are close to one. Thus the error introduced by 
this simplification is small relative to the precision 
of our computation. 



• MINOS: The MINOS experiment utilizes 120 
GeV protons from the NuMI beamline impact- 
ing a graphite target. The near-detector has a 
large overall mass, but only part of the cross 
sectional area of the detector is instrumented 
and the near-detector itself is nearly 1 km from 
the target. Nonetheless, the large boost pro- 
vided by the 120 GeV proton beam leads to 
a significant event rate for dark matter scat- 
tering. We use an estimate for the total POT 
prior to the 2012 shutdown for the NOvA up- 
grade. 

• T2K: The T2K experiment has only been op- 
erating since 2010, and rather than use the 
current dataset we have taken into account 
the final number of POT expected for the 
current run. T2K utilizes a 30 GeV proton 
beam impacting a graphite target, and has 
two near detectors, ND280 and INGRID, both 
located in a complex 280 m from the tar- 
get. ND280 is about 2 degrees off- axis and is 
well-instrumented with TPCs for tracking and 
analysis, while INGRID is on-axis and contains 
significantly more fiducial mass. 

• MiniBooNE: The MiniBooNE experiment uti- 
lizes an 8.9 GeV proton beam impacting a Ger- 
manium target and, distinct from MINOS and 
T2K, has a single spherical mineral oil detec- 
tor located 541 m from the target. This de- 
tector has a large fiducial mass, and impor- 
tantly MiniBooNE has already published a full 
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FIG. 7. Tree-level dark matter scattering off nucleons mediated by the vector and scalar portals. 



neutrino elastic scattering analysis [19] with 
O(IO^) events and a measured energy spec- 
trum, which provides the natural background 
for any dark matter beam search. We use an 
estimate for the total POT prior to the 2012 
shutdown for the NOvA upgrade. 

The simulation of the dark matter beam used a 
re- weighting technique, first determining the dark 
matter trajectories that intersect the detector, and 
subsequently weighting them according to the pro- 
duction distributions discussed earlier in Sec. [21 We 
will describe these two steps in more detail below, 
starting with the generation of the dark matter tra- 
jectories. 

For direct production at either T2K or MINOS, 
the y's were generated over an array of kinemat- 
ically allowed momenta, and each V was decayed 
isotropically into a random pair of x's in the V's 
center of mass frame. The lifetime of the V is short 
enough for the parameter space considered that it 
will decay before escaping the target, and so the 
propagation of the V through the target is ignored 
in the simulation. The trajectories of each of the x 
particles are then checked to determine if they pass 
through the fiducial volume of the corresponding 
near detector. These trajectories are recorded along 
with the energy of the x- The treatment of indirect 
production at T2K, MINOS and MiniBooNE was 
similar (see [H]), but required the extra initial step 
of first generating kinematically allowed meson tra- 
jectories, with each then decayed isotropically into 
a V and a 7 in the meson rest frame. The newly 
produced V is then treated in the same manner as 
in the direct production simulation. 

With the trajectories in hand, for each point in pa- 
rameter space the expected number of events could 
be determined by weighting them according to the 
production distribution /(6>,p), the scattering cross 
section aJ^^{E), and the distance R which x propa- 
gates through the detector. There is also an overall 
measure factor: A = 5p50S(j)/{27T) for indirect pro- 
duction, OT A = dp for direct production, where the 
6 quantities refer to the step sizes used in the sim- 
ulation for (j) or V production. Note that the dis- 
tance R travelled through the MINOS near detector 
and ND280 will almost always equal the length of 
the detector Ldet shown in Table |I| For INGRID, 



it will occasionally be twice the listed number if it 
passes through the center of the detector, where two 
of the detector's modules overlap. MiniBooNE uses 
a spherical detector, and so R can vary significantly 
in this case. 

The final expression for the expected number of 
elastic nucleon dark matter scattering events is given 

by 



Nx^Nx = riN X Ceff 



(25) 



E 

prod, 
chans. 



E 

trajec. 



R,af^{E,)f{0,,p,)A, 



where tiat is the nucleon density in the detector, 
while eeff is the detection efficiency for events within 
the specified fiducial volume and cuts on momentum 
transfer. We will assume that lower cuts are above 
the range for coherent elastic scattering, so that our 
nucleon-level treatment in (24) should be reliable. 



We will also assume that the detection efficiencies 
do not deteriorate significantly for the full range of 
momentum transfer relevant for DM scattering. The 
production quantities are given by 



2A^POT ^ utIt^ pt 
2N^ X Br((/:) ^ X + 



direct 
indirect 



(?6) 



fv{p) X 1(1 
/™^(^,P) 



COS' 



9) 



direct 
indirect 



(27) 



The distributions for direct {fv{p)) and indirect 
{f^^{0,p) = f^^^^{0,p) or f^iO.p)) production 



were discussed in Sec. 2|^ Note that the meaning of 



p and 6 varies depending on the context. For direct 
production, p is the V momentum, and is the an- 
gle between the dark matter and the beam in the 
V rest frame. For indirect production, both p and 
refer to those of the original meson in the lab 
frame. The direct production parameters in are 



2 For T2K, rather than f^^^^(0,p), the indirect production 
distribution used was a parametrization of data from NA61 
20 , using a replica T2K target. However, the results are 
consistent with those using the BMPT parametrization. 
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Target 


It 


POT 


-E'beam 


L 


Adet(cm^) 


Ldet 


nAr(cm ^) 


Fiducial Mass 


Ceff 


MiniBooNE 


Be 


71 cm 


1.2 X 10^^ 


8.9 GeV 


541 m 


1.2 X 10^ 


11.5 m 


9x10^^ 


^650 tons 


0.6 


MINOS 


C 


94 cm 


1.5 X 10^^ 


120 GeV 


965 m 


7.1 X 10^ 


1.3 m 


5x10^^ 


27 tons 


0.8 


T2K ND280 


C 


90 cm 


5 X 10^^ 


30 GeV 


280 m 


5.5 X 10^ 


0.7 m 


4x10^^ 


1.67 tons 




T2K INGRID 


C 


90 cm 


5 X 10^^ 


30 GeV 


280 m 


2.2 X 10^ 


0.585 m 


5x10^^ 


--110 tons 





TABLE I. A summary of the parameters used for the three experiments considered in this work; see e.g. [191 I35ti37| for 
MiniBooNE, [38ti42] for MINOS, and [43^46] for T2K. Further details are in the text including a description of the notation. 
Note that in the absence of published analyses focusing on neutral currents, the overall efficiency Ceff for T2K is not known; we 
take it to be of the same order as Ceff for MINOS. 



the number of protons on target A^poT, the target 
length It and density tit, and the total cross section 
cTpT- The experimental quantities are listed in Ta- 
ble |l| while our treatment of apx was discussed in 
Sec. [2] For indirect production, we use an estimate 
of the total (j) = rj yield Nr^ and the branching ratio 
to the mediator. 



C. Results 

The results for DM nucleon scattering at MINOS, 
T2K, and MiniBooNE are shown for various parame- 
ter choices in Figs.[8][l3j All the plots show contours 
of the number of events (10, 1000, or 10^) in the 
plane of nucleon scattering cross section (or kinetic 
mixing n^) versus dark matter mass. The sensitiv- 
ity tends to be fairly flat as a function of m^, as 
the momentum transfer in the scattering tends to 
be much larger than the mass and thus drops 
out of the kinematics. The exception to this general 
rule is that when approaches the decay thresh- 
old, ~ my/2, there is an enhancement in sen- 
sitivity as the dark matter has a small transverse 
boost from the V decay and thus a larger fraction 
of trajectories will intersect on-axis detectors. The 
overlayed dark line denotes the parameter choices 
consistent with x having the correct thermal relic 
density to form WIMP dark matter. The structure 
in this curve for higher mass reflects the p/uo and (j) 
hadronic resonances that play a role in annihilation. 
Note that the position of the thresholds and reso- 
nances is shifted down slightly from 2m^ due to the 
kinetic energy of the WIMPs. Assuming an initial 
thermal abundance, the WIMP relic density would 
be too large in regions of parameter space below this 
curve. 

Existing particle physics limits on the parame- 
ter space would, as discussed in Sec. [3) require 
that n and a' satisfy both {a/a')^!^ < 10~^ and 
a' < few X 10~^ for models of this type. For 
the perturbative values of a' ^ a that will pri- 
marily be used in what follows, the former limit 
is more restrictive implying K:'^ < 10~^. However, 



this constraint is inferred from 5-factory analyses 
of somewhat different models rather than dedicated 
searches, so we refrain from showing any explicit ex- 
clusion curves. Nonetheless, with this benchmark 
in mind we observe from the plots that interesting 
sensitivity emerges with the ability to distinguish 
O{10^ — 10^) dark matter scattering events from the 
neutrino background. The characteristic scattering 
cross section per nucleon to achieve (9(1000) events 
ranges from 1 — 10 pb, which is an impressive level of 
sensitivity. Since this is not coherent scattering, it 
is more useful to contrast it with the best low-mass 
sensitivity achieved for spin-dependent scattering in 
underground detectors [47] which is around O.lpb for 
^ 10 GeV, but drops off rapidly for lower masses. 
Of course, to implement a search for (9(1000) events 
would require the ability to separate this from the 
neutrino background of (9(10^ — 10^) events. We will 
comment further on possible search strategies and 
means for background rejection in the next section. 

The individual plots reveal a number of other fea- 
tures summarized below: 

• In Figs. [S] and [9] we exhibit the sensitivity to 
dark matter in the 100-200 MeV mass range, 
with mv = 400 MeV. The direct parton-level 
approximation for on-shell production of such 
low mass vectors is questionable with the PDF 
scale set to Q = my, so we use indirect pro- 
duction via 7^-decays. As there is no signif- 
icant resonant enhancement in r] production, 
one should bear in mind that this may be an 
under-estimate for total F-production in this 
mass range. The ensuing sensitivity is shown 
for T2K in Fig.^ and for MiniBooNE and MI- 
NOS in Fig.[9j Note that the lower momentum 
transfer in scattering at MiniBooNE leads to 
an enhanced sensitivity to the cross section for 
a given sensitivity in relative to MINOS. 



• In Figs. [10] and [IT] for T2K and Fig. [12] for 

MINOS we exhibit the sensitivity to higher 
mass dark matter, for my = 1 and 2 GeV 
using direct parton-level production. We ex- 
pect the narrow width approximation to work 
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FIG. 8. Expected number of neutral current-like dark matter nucleon scattering events from V's produced through 77 decays 
for the ND280 (left) and INGRID (right) detectors at T2K with my = 400 MeV. The regions show greater than 10 (hght), 1000 
(medium) and 10^ (dark) expected events. The dashed curve indicates the value of k required for the dark matter annihilation 
cross section in the early universe to equal 1 pb. 
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FIG. 9. Expected number of neutral current-like dark matter nucleon scattering events from V's produced through 77 decays 
for the MINOS near detector (left) and MiniBooNE (right) with my = 400 MeV. The contours are described in Fig. [s] 



fairly well for 2 GeV vectors with this mass set- 
ting the PDF scale. The precision of the esti- 
mate will certainly be lower using this method 
for 1 GeV vectors, due to the uncertainties in 
the PDFs and the importance of higher-order 
QCD corrections. Nonetheless, we see that the 
tree-level sensitivity for 1 GeV vectors is only 
marginally enhanced relative to my = 2 GeV. 



Using Figs. 10 and[TT] it is interesting to com- 
pare the sensitivity of the two near-detectors 
at T2K. Given the suppression of direct X" 
production in the forward direction, the off- 



axis ND280 detector at T2K is ideally posi- 
tioned to capture a comparatively large flux of 
dark matter, as compared to the on- axis de- 
tector INGRID. However, the much larger ac- 
tive mass of INGRID more than counteracts 
this effect, leading to an enhanced sensitivity. 
We are not aware if T2K has plans to use IN- 
GRID for analyses unrelated to diagnostics of 
the neutrino beam, but we see that there is 
considerable intrinsic sensitivity to light dark 
matter. 



• In Fig. 13 we compare the direct production 
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mv = 1 GeV a' = a mv = 1 GeV = a 
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FIG. 10. Expected number of neutral current-like dark matter nucleon scattering events from direct V production for the 
ND280 (left) and INGRID (right) detectors at T2K with my = 1 GeV. The contours are described in Fig. [s] 
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FIG. 11. Expected number of neutral current-like dark matter nucleon scattering events from direct V production for the 
ND280 (left) and INGRID (right) detectors at T2K with my = 2 GeV. The contours are described in Fig. [s] 



sensitivity of the INGRID detector at T2K for 
two values of the dark U(l) couphng = a 
and a' = 1. The latter value implies a self- 
interaction cross section for dark matter a ~ 
47rm^/my that can reach (9(0.1 mb) for a 1 
GeV mediator. This is close to, but somewhat 
below, the scale that would lead to detectable 
effects on halo structure, which may be rele- 
vant to the understanding of the inner regions 
of dwarf spheroidal halos. As is apparent from 
the plot, increasing has the effect of enhanc- 
ing the annihilation rate and thus moving the 
relic density curve to lower values of hi. Conse- 
quently, this increases the intrinsic sensitivity 
to K while effectively lowering the sensitivity 



to the scattering cross section. Analogous sen- 
sitivity to applies to the other parameter 
regimes and experiments shown in the earlier 
plots which all assume a' = a. 

• The plots all indicate that very light WIMPs 
with masses below about 100 MeV are prob- 
lematic as thermal relics. Even for larger val- 
ues of a', the event rate along the measured 
relic density curve grows as decreases, and 
reaches levels which are well above the elastic 
scattering rate for neutrinos. Thus, the mea- 
sured elastic scattering of neutrinos at these 
facilities, and its consistency with the Stan- 
dard Model, serves to exclude a large class of 
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FIG. 12. Expected number of neutral current-like dark matter nucleon scattering events through direct V production for the 
MINOS near detector with two different vector mediator masses (mv = 1 GeV on the left and mv = 2 GeV on the right). The 
contours are described in Fig. [s] 
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FIG. 13. Expected number of neutral current-like dark matter nucleon scattering events from direct V production with the 
INGRID detector at T2K, comparing two different a values (o;^ = o; on the left and o;^ = 1 on the right) for a 2 GeV Vector 
mediator. The contours are described in Fig. [s] 



models of MeV scale dark matter. This strong 
tension with models of MeV-scale dark mat- 
ter was already exhibited in more detail using 
data from LSND and MiniBooNE in 



5. CONCLUDING REMARKS 

The direct search for dark matter is above all 
a search for weakly-interacting degrees of freedom, 
and at this stage the mass range is relatively uncon- 
strained. The LHC has as yet revealed little sign of 
new weak-scale physics, so it is important to keep 



in mind that the simple thermal relic paradigm is 
broad enough to encompass a large mass range, ex- 
tending well below the weak scale. Thus it is crucial 
to utilize all the available experimental tools to ex- 
plore the viable dark matter parameter space. The 
weak nature of DM-SM interactions means that fixed 
target neutrino experiments provide a very natu- 
ral source of low mass sensitivity. This goes both 
ways, as the next generation of underground dark 
matter direct detection experiments may in turn be 
able to detect various astrophysical and cosmolog- 
ical sources of neutrinos. Rather than being sim- 
ply an irreducible background, it seems clear that 
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these experiments will have to become observatories 
for all types of cosmic weakly-interacting degrees of 
freedom. In this paper, we have explored another as- 
pect of this convergence, namely the use of neutrino 
beam experiments to probe light dark matter that 
can be produced in the target, and undergo elastic 
scattering in the detector. 

The challenge in developing a search strategy for 
this dark matter signal will be in disentangling the 
event spectrum from the neutrino background at 
a level of maybe 1-10%. This would allow sensi- 
tivity to kinetic mixing in the f^^ < 10-"^ - 10-^ 
range, which is the most viable regime given the 
level of indirect constraints. This level of sensitiv- 
ity to sub-GeV WIMPs, at the pb-level in terms 
of per-nucleon scattering cross section, is only at- 
tained in spin-dependent direct detection for WIMPs 
with much larger masses, exceeding 10 GeV. Thus 
neutrino experiments could provide an important 
means of probing dark matter nucleon scattering be- 
low the mass range accessible via direct detection. 
In terms of isolating dark matter beam scattering 
events from the large background of neutrino elastic 
scattering, we note that there are several distinc- 
tive characteristics. Firstly, the dark matter beam 
has a higher average energy (shown in Fig. [5| than 
the neutrino beam, and in particular a much higher 
cutoff that approaches the energy of the primary 
proton beam. This would permit a relatively high 
cut in momentum transfer in scattering, provided 
such events are retained in the full sample. Sec- 
ondly, the dark matter beam will be relatively unaf- 
fected by turning off or switching the polarity of the 
magnetic focusing horns, which would alter the neu- 
trino beam significantly. Finally, there may also be 
useful information in the (nanosecond-scale) timing 
structure as the production mechanisms for vector- 
portal-coupled DM and the neutrino beam are quite 
distinct. Determining whether one or more of these 
features could be put to practical use in a search 
strategy would require a dedicated analysis. 

In concluding, we would also like to comment 
on some alternative approaches to explore the light 
WIMP regime. 



• Direct detection in the low mass range could be 
feasible using electron scattering, as explored 
in recent work [48]. This approach is quite 
complementary to the neutrino beam analy- 
sis considered here. While the beam analysis 
requires relatively heavy vectors with my > 
2m^, the sensitivity for electron scattering is 
enhanced when my <C m^. For comparison, 
the projected electron scattering sensitivity to 
the vector portal model considered here is rela- 
tively weak for my ~ 1 GeV, but becomes sig- 
nificant for my ~ 1 MeV [48] . Future progress 
using Ge crystals seems promising and may al- 
low strong sensitivity to sub- GeV WIMPs with 
very light MeV-scale mediators, provided tech- 
niques are available to deal with all the back- 
grounds at such low recoil energies. 

• Direct collider searches are also possible, 
utilizing missing energy signatures such as 
monophotons or monojets [49 . However, the 
sensitivity weakens significantly with light me- 
diators. Collider searches can also pursue sig- 
natures of the light mediators directly, e.g. via 
subleading SM decays which may produce sig- 
natures such as lepton jets at high energy. In 
the scenarios considered here, the decays of the 
mediator are all prompt so there are no dis- 
placed vertices. 

Looking to the future, the continued development 
of long-baseline neutrino facilities provides an ideal 
setting for expanding the search for light hidden 
sector states. This sensitivity extends beyond the 
models of light dark matter discussed here to other 
classes of new physics, such as the scenarios dis- 
cussed in [50] , that could impact the neutrino sector 
more directly. 
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